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ABSTRACT 

We consider the pinning of superfluid (neutron) vortices to magnetic fluxtubes associated with a 
type II (proton) superconductor in neutron star cores. We demonstrate that core pinning affects the 
spin-down of the system significantly, and discuss implications for regular radio pulsars and magnetars. 
We find that magnetars are likely to be in the pinning regime, while most radio pulsars are not. This 
suggests that the currently inferred magnetic field for magnetars may be overestimated. We also 
obtain a new timescale for the magnetic field evolution which could be associated with the observed 
activity in magnetars, provided that the field has a strong toroidal component. 
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1. INTRODUCTION 

The outer core of a mature neutron star is expected 
to contain a neutro n superfluid couple d to a type II pro- 
ton superconductor (jBavm et al.lll969f ). The dynamics of 
this system is complicated; the neutrons rotate by form- 
ing a dense array of vortices while the magnetic field 
is carried by quantized proton fluxtubes. The interac- 
tion between the vortices and the fluxtubes is expected 
to have significant impact on the dynamics of these sys- 
tems. In fact, the vortices may "pin" to the fluxtubes 
(in the sense that the energy cost of cutting through is 
too high) as the star spins down, affecting the evolution 
of the star's rotation and magnet ic field. An example of 
this was provided by iLinkl (|2003l ). who argued that the 
free precession interpretation of the long-term variabil- 
ity in PSR BI828-II would not be compatible with core 
vortex pinning. Vortex-fluxtube pinning may also cou- 
ple the magnetic field evolution t o the rotation, as in the 
model of iRuderman et ail (|1998l ). 

Connections between the magnetic field and the spin 
evolution are crucial if we want to understand the dif- 
ferent observational manifestations of neutron stars, like 
radio pulsars and magnetars, and establish an evolution- 
ary link between them. We also need to keep in mind 
that the spin and its rate of change are the primary ob- 
servables. The magnetic field tends to be inferred from 
the assumption of pure dipole breaking. 

The aim of this Letter is to demonstrate how the stan- 
dard logic linking the spin evolution and the inferred 
magnetic field breaks down if there is significant vortex- 
fluxtube pinning, and how the same mechanism could 
drive the magnetic field evolution. 

2. A VORTEX DYNAMICS MODEL 

We consider the motion of individual neutron vor- 
tices and proton fluxtubes in a neutron star core, com- 
posed of neutrons, protons and electrons. The local vor- 
tex/fluxtube velocity is u\ where the index x = {n, p} 
distinguishes between the two types. For simplicity, the 
hydrodynamical fluid velocities are assumed to represent 
rigid-body rotation, i.e. t> x = e lJ CtjXk, where the spin 



frequencies J7 X are aligned (defining the z-axis) but may 
differ in magnitude. A direct consequence of uniform 
rotation is that the neutron vortices are straight, which 
may be a serious oversimplification. The proton flux- 
tubes are not constrained in this way. 

In order to make progress, we make a number of simpli- 
fications. The most basic is the assumption of axisym- 
metry. We ignore the entrainment between the super- 
fluids on the hydrodynamic scale [although the effect is 
inclu ded implicitly, in th e magnetization of neutron vor- 
tices (jAlpar et alJ ll984)]. Given these assumptions one 
can use Ampere's law, as modified in a type II super- 
conductor, to show that the electron and proton fluids 
are "locked" together, movin g with a common velocity 
Vp (jGlampedakis et all l2011h . Finally, we identify fl p 
with the motion of the "crust" and the observed spin 
frequency Q. 

The motion of an individual vortex/fluxtube is gov- 
erned by a force balance equation. Ignoring the inertia 
of the filaments (as usual) we have 



fun + /dh + /pn — 
/mp + /d p + fp P + /t — 



(1) 
(2) 



Each vortex/fluxtube experiences (i) a Magnus force 
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p x e ijk K]{u" k 



'£) due to the interaction between 
the macroscopic fluid flows and the mesoscopic, quan- 
tized, neutron/proton circulations (p x is the fluid density 
while kJ; has magnitude k = h/2m and points along the 
local vortex/fluxtube direction), (ii) a drag force = 
p x nlZ x (Vp — u x ), typically attributed to the scattering of 
electrons by the intrinsic vortex/fluxtube m agnetic field, 
with a dimensionless coeffic ient 1Z X <C 1 (Alpa r~e"t al.l 
11984 IRuderman et al.|[l998fl and (iii) a "pinning" force 
/p x due to local interaction between a vortex and a flux- 
tube. Finally, a bent fluxtube will experience a (self- 
induced) tension [since we are considering straight 
vortices there is no similar force in (JTJ]. All forces are 
per unit length. 

The direct interaction is primarily due to the 
short range magnetic fields associated with each vor- 
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tex/fluxtube, and could lead to pinning dSauli 1X9891. 
This interaction must obey Newton's third law, i.e. 



/] 



p P 



-«/pn 



(3) 



In fact, this relation is the operational definition of /p . 
In a typical fluid element the number density (per unit 
area) j\f p of fluxtubes greatly exceeds the vortex density 
A/" n , resulting in only a tiny fraction of fluxtubes inter- 
acting with vortices at any given instant. At the same 
time neighbouring fluxtubes can push each other, "dis- 
tributing" the action of the vortex back-reaction force 
— /p n to the entire fluid element. On average, this effect 
can be quantified by the small paramet er a ~ j\F n /j\f p 
(|Ruderman et alJll998t[Jahan-Mirill2000t l. 
Eqn. (JT|) provides an exact form for the pinning force, 

/ Pn = -p n nU n (vl - <) + p^ k z, {vl - ul) (4) 

where we have set n l n = nz l . This can be used in ([2]) to 
eliminate / Pp [with the help of eqn. ([3])]. Thus we obtain 

TZ(v; - <) + e^ Zj (ul - V k ) +T=0 (5) 
where we have defined 



Up 
K 
1 

= K 



V 1 = 



K\\e ijk Zj uf l 



(6) 
(7) 



with /C = km + a/xp. We have introduced the rela- 
tive vortex- fluxtube velocity w* cl = u l n — ut, the ratio 
x p = Pp/ p-a (essentially the proton fraction) and the (di- 
mensionless) projections km and k± such that 



KHZ + K | 



1 



(8) 



Given that the macroscopic magnetic field B l is locally 
directed along kt, the decomposition ([8]) makes contact 
with a general "twisted-torus" magnetic field configu- 
ration which has a mixed toroidal (-Bt) & n d poloidal 
(Bp) character represent ative of hydromagnetic equilib- 
rium (jBraithwaitel 12009V Then K|| ~ B P /B T and for 
a typical twisted-torus configuration km ~ k± ~ 1. A 

predominantly toroidal field has km <C 1 and n\ ~ tp l . 
Excluding the special case km < oi/x v of a (essentially) 
purely toroidal field, we assume K. ~ km hereafter. 

As a final step, we can invert ([5]) and obtain the vortex 
velocity 
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3. THE PINNING REGIME 



Strong interaction between a vortex and a fluxtube seg- 
ment (via the forces /p x ) may result in more or less per- 
fect pinning. Given the present understa nding of neu - 
tron star matter, this is likely to happen (|Saulsl[X989D . 



Assuming vortex pinning, in the remainder of this Let- 
ter we explore some very interesting consequences. If 
(for whatever reason) pinning is inefficient, in the sense 
of a significant vortex-fluxtube relative motion, then our 
model does not apply. 
When pinned, the two segments would move together, 

< — Up = (10) 

This velocity can be immediately obtained from ([9]) after 
setting u l rcl = 0. Using cylindrical coordinates {w, ip, z} 
with respect to the rotation/symmetry axis, we find 
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(ii) 



where we have defined the rotational lag Q np = O n — £1. 
We see that, to leading order, the velocity of pinned vor- 
tcx/fluxtube segments is nearly azimuthal and in corota- 
tion with the charged component. Superimposed on this 
there is a small relative motion Au l = u 1 — v p between 
the vortex/fluxtube array and the proton velocity with 
both radial and azimuthal components: 



K|| 
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Kt 



fro 
JT 
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Au v 



(12) 
(13) 



The dependence of Au l on Ky encodes the action of the 
(radial) Magnus force on a fluxtube. We note that Ait 1 
is nearly azimuthal unless the magnetic field is dom- 
inated by the toroidal component to the degree that 



io- 



The true magnitude of 1Z P is un- 



certain. We use a valu e inferred from ele ctron scattering 
by a sin gle fluxtube (l Alpar et al.lll984h . However, we 
note that Uonesl (|2006f ) has argued that the high fluxtube 
density suppresses the magnetic scattering leading to a 
much smaller 1Z P . Nevertheless, as no actual alternative 
value for 1Z P has been provided in the literature, we base 
our discussion on the standard "single fluxtube" result. 

A vortex cannot move arbitrarily fast and remain 
pinned. Above some relative vortex-fluid velocity thresh- 
old the forces entering the balance ((H) will exceed the 
maximum pinning force, / p i n , and the vortex unpins. 
The maximum p inning force has been estimated to be 
[e.g.EEajlQQl], 

/ pin pa 3 x 10 15 Bli 2 dyn/cm (14) 

where B\2 = -B/10 12 G is the normalized interior mag- 
netic field. At maximum pinning, the vortex force bal- 
ance is well approximated by after setting u n « u l p « v p 
and omitting the drag term in (j4|)) 



(15) 



and we can extract the maximum allowed spin lag (above 
which pinning can no longer be sustained); 



q m ax 
5 : np 



/pin 
p n KU7 



7.6 x 10" 



B 



1/2 



TU 6 



(16) 



where zuq = G7/10 6 cm. This result can be used in 
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(|12p and (|13|) to provide the maximum relative vor- 
tex/fluxtube velocity with respect to the protons; 



(a/pi„ + /?) (17) 



Let us provide numerical estimates of the forces in (I17p . 
The tension can be approximated by /t = E^jR c where 
is the fluxtube energy per unit length and R c is the 
lo cal radiu s of cu rvature. Estimates for can be found 
in Mcndclj (|1991[) while a typical value for the curvature 
radius would be R c ~ R. The pinning force can be esti- 
mated using (|T4| together with a = Af n /Af p , 7V n = 2VL n / n 
and A/p = B/<po (<fo = hc/2e is the flux quantum). This 
way we obtain 



10~ 9 cm/s, 



10 



-10 



PB 



V2 n 



- cm/s (18) 



with the stellar spin period, P, measured in seconds. 

4. SPIN EVOLUTION 

Vortex pinning has direct impact on the neutron star 
spin evolution. This is not surprising, given the close link 
between the spindown of the neutron superfluid and the 
radial vortex motion, 



-CC7il n /2f2 n 



(19) 



If the neutron superfluid is able to track the 
electromagnetically-driven spindown of the charged com- 
ponent (i.e. when 57 n = D, and f2 n = f2), we have 



8 x lQ- 7 m 6 (I0 4 yr/r sd ) cm/s 



(20) 



where r s( j = P/2\P\ is the observed "spin-down age". 
This vastly exceeds the maximum radial vortex velocity 
(|17[) in the pinning regime. Thus, we should expect the 
neutrons to spin down at a rate | f2 n | <C | f2 in the pinning 
regime. This would lead to a monotonically increasing 
spin lag f2 np , to the point where f2]^f x is reached and 
vortex pinning can no longer be sustained. The develop- 
ment of this lag is similar in most models of large pulsar 
glitches. 

Let us focus on the timescale associated with the pin- 
ning phase. In order to estimate this, we need to model 
the spin evolution of a multi-fluid star. This can be done 
by considering the volume-integra ted Euler equations for 
the fluids [cf. iSiderv et all (|2010t )]. The resulting equa- 
tions take the form, 



I n fl n = N m f, 



I p Sl = -N m{ 



(21) 



where J x = J p x zu 2 dV is the moment of inertia for 
each fluid. These equations feature torques due to (i) 
the exterior magnetic field (iV* m ) and (ii) the "mutual 
friction" (iV^ f ) which represents the vortex/fiuxtube- 
mediated coupling. The latter torque is given by 



N, 



mf 



e ijk Xj F^ f dV 



(22) 



where the mutual friction force is 



*£f = A/"n (fhn + fk) = 2n nPn e^ k z, (vl - u k ) (23) 



Using our result (|lip for the common vortex/fluxtube 
velocity in the pinning regime we have 



Nli « -2Q a z' l w 



II 



Tt 

VjppK 



dV (24) 



We also need to consi der the magnetic t orque 1\ 
Following the analysis of [Spitkovskyj (|2006l ) , which 
counts for a wind contribution, we have 



cm ' 

ac- 



N l = 

cm 



gig 

4c 3 



n 3 z l = Ati 3 z l 



(25) 



where P>d is the exterior dipole field at the magnetic pole. 

Returning to the spin evolution equations (|21|) with 
(|24|) and p5|) . we observe that the timescale r m f ~ 
I n fl/N m f associated with the pinning and tension forces 
is much longer than the electromagnetic spindown 
timescale. For typical neutron star parameters we find 
Tmf ~ 10 10 k| yr. Hence, the spin evolution in the perfect 
pinning regime is well approximated by 



Cl n W 0, Q « -Mf/Ip 



(26) 



in accordance with the intuitive notion of a limited radial 
vortex motion enforced by pinning, and of the magnetic 
braking affecting only the charged component. Solving 
(f26|) with initial data O x (0) is straightforward, and we 
have 



n(t) = fi(0) 1 



Ip TO 



-1/2 



fin(0) (27) 



where Iq = I n + I p is the total moment of inertia and 
To = Iq/(2A£1 2 ). In a single-component star, tq would 
correspond to the spin-down timescale and would be 
identified with the observed spin-down age r s d. How- 
ever, in the superfluid pinning model the appropriate 
identification, based on the spindown law (|26|) . is r s d = 
(I p /Io)tq sa x p tq. Equivalently, we can relate the true 
field R>d with the field B d inferred from spin-down, 



1/2 



7T 2 i? 6 



\P\P 



1/2 



r l/2 R inf 
x p n d 



(28) 



The spin evolution (|27[) naturally leads to a monotoni- 
cally increasing spin-lag. Considering, for instance, the 
early regime t -C (I p /Iq)tq, we can easily obtain the time 
t max at which the system reaches the maximum lag (|16p , 
where vortices first unpin [assuming O np (0) = 0] 



tr 



oQmax 
zl 'np 

n 



2.4 x 10 



pol/2 
-3 - r£> 12 

zu 6 



(29) 



which features the interior B field. When applied to a 
typical radio pulsar (P = 0.1s, B = 10 12 G) the result 
|29l) leads to t 

max <K t sc [. Thus, core vortex pinning in a 
typical pulsar is an ephemera l phenomenon, plausibly as- 
socia ted with large glitches (jGTampcdak is fc Anderssonl 
l2009h . 

The situation may be very different for magnetars. 
For a strongly magnetized and slowly spinning object 
the timescale t m ax can be comparable to, or even ex- 
ceed, r s d- In this case we need to resort to the full 
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solution (|27l) . Inserting canonical magnetar parameters 
(B = 10 15 G, P = 10 s), we obtain 

imax ^ 7"sd (30) 

This suggests that vortex-fluxtube pinning may well per- 
sist in magnetar cores. 

5. IMPLICATIONS AND DISCUSSION 

We have provided a rather simplistic analysis of a 
very complex problem, Yet, the results could have im- 
portant implications for the magneto-rotational prop- 
erties of neutron stars. Perhaps the most significant 
conclusion concerns the potential persistence of vortex 
pinning in magnetars, as suggested by (|30"|) . This im- 
plies (see eqn. ([25)) ) that the surface dipole field in- 
ferred from spin-down is an overestimate by a factor 

(I-p/Iq)^ 1 / 2 ~ Xp 1 ^ 2 w 3 — 10. This mechanism could 
be responsible for the apparent clustering of the ob- 
served magnetar spin periods, since the end of the pin- 
ning regime would be accompanied by an effective weak- 
ening of the torque; N cm — s> x p N em . Interestingly, the 
presence of a "twisted " magnetar magnetosp here may 
lead to a similar effect ([Thompson et alll2002l ). 

There are also implications for the magnetic field evolu- 
tion in magnetar cores. According to the pinning model 
the motion of the proton fluxtubes (and therefore of the 
magnetic field) is controlled by the neutron vortices and 
vice versa. We can define a magnetic evolution timescale 
associated with the relative vortex/fluxtube motion with 
respect to the proton fluid, tb = £/|Ait max |, where L is 
some typical lengthscale. In magnetars Aw max is mostly 
due to the tension since (fT8|) implies /t S> a/ P i n [note 
that (Q~6|) suggests il n ~ 57 for magnetars]. The short- 
est timescale is the one associated with the azimuthal 
motion Ait^ ax , leading to 

K\\np p L 6 Bp f R c \ 
tb~^—~10 Li—^—jyr (31) 

where we have restored the explicit dependence of /t on 
the curvature radius R c and defined L5 = L/10 5 cm. The 
timescale associated with the radial fluxtube motion is 
longer by a factor 1Z~ 1 ~ 10 3 . Also, tb becomes signifi- 
cantly longer if (for whatever reason) /t is negligible. 

Is the timescale in (f3lT) relevant for the observed mag- 
netar population? For this to be the case we need tb to 
be comparable to the age r s( j ~ 10 3 — 10 4 yr of these ob- 
jects. At first sight, (I5T|) seems to suggest a significantly 
longer timescale. However, our result approaches the in- 
teresting regime if we consider a magnetic field with a 
locally strong toroidal component Bt > -Bp and a cur- 
vature R c < O.li?. In deed, recent work on stable hy- 
dromagnetic equilibria (Braithwaitc 2009) suggests that 
these properties could be realistic. 

Our results indicate that in the regime of persistent 
vortex pinning in magnetar cores the evolution of a 
strong toroidal magnetic field component is accelerated 
with respect to that of the poloidal component, and could 
occur on a timescale 10 4 — 10 5 yr. Comparing to recent 
work on ambipolar diffus ion in superfluid neutron stars 
(jGlampedakis et al.|[20TTI ). we see that persistent vortex 
pinning is the main driving force of magnetic field evo- 
lution in magnetar cores, with possible relevance to the 
activity observed in these objects. We also note that tb 



can be comparable to the ambipolar diffusi on timescale 
in we akly / non-superfiuid magnetars (e.g. lArras et al.l 
(2004)). ' 

The implications for radio pulsars are more subtle. Re- 
calling equation (|29[) . core vortex pinning may only be 
supported over short periods of a pulsar's life. Thus, 
pulsars likely spend most of their time in a "cutting 
regime" wh ere vortices can move through the fluxtubes 
(lLinkl l2003D. This result is in ap pare nt conflict with th e 
work bv lRuderman et al.l (|1998f) and I Jahan-Mhl (|2000lh 
In these studies it is assumed that the radial vortex mo- 
tion is given by (|19p with Q a = Q and Sl n = f2 in the 
pinning regime. Our discussion shows that this assump- 
tion is not valid. Moreover, the previous studies do not 
take full account of the Magnus force. Hence , the results 
of IRuderman et al.l (|1998D and IJahan-Minl (|2000l) must 
be considered with caution. 

In this Letter we have discussed issues that are central 
to our understanding of the magneto-rotational evolu- 
tion of neutron stars. We have shown that the observed 
magnetars may be in the regime of core vortex pinning, 
affecting the bulk spin-down, the inferred surface mag- 
netic field and the magnetic field evolution. Our results 
also show that core pinning is likely less relevant in the 
weaker magnetic field radio pulsars. 

The issues we have discussed obviously need more de- 
tailed consideration. Future work needs to relax many 
of our assumptions. First of all, we have assumed that 
the entire star is in a superfluid state while, in real- 
ity, this may be an accurate description only for the 
outer core - in that case Bd would be overestimated 
by the smaller factor ~ ((/ p + In)/Io) where In is 
the moment of inertia of the non-superfiuid neutrons. 
To make progress we need to consider a model with re- 
alistic stratification. For the magnetars in particular, 
it not not yet clear that their core temperature is be- 
low the level for the onset of neutron superfluidity sug- 
gested by the recent coolin g mod els for the Cas siopei a 
A neutron star ()Page et all 120 lit iShternin et~aTll2011l) . 
The most detailed work on magnetar temperature pro- 
files (jKaminker et al.l r2006f) appears t o favour superflu- 
idity, but the issue is far from settled (jArras et al.1120041 ; 
IDall'Osso et~a l. 2009). A related case for which our 
model does not apply is for a bulk B field in excess of the 
critical value ~ 10 16 G above which pr oton superconduc- 
tivity is suppressed ([Bavm et al.lll969l ). We also need to 
remove the assumption of uniform rotation, which does 
not remain consistent when the long-term neutron star 
evolution is considered. This would help link our vor- 
tex/fluxtube model to the physics of ambipolar diffusion 
(Gla mpedakis et aL1l2011|) . The departure from uniform 
rotation is also likely to be important in the vortex cut- 
ting regime, which should be relevant for normal pulsars. 
These issues require us to consider more realistic neutron 
star models, a challenge that should ultimately reward us 
with a deeper understanding of these exciting systems. 
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